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We investigate the effect of a vertical electric field on a Dirac semimetal thin film. We show 
that through the interplay between the quantum confinement effect and the field-induced coupling 
between sub-bands, the sub-band gap can be tuned and inverted, during which the system undergoes 
a topological phase transition between a trivial band insulator and a quantum spin Hall insulator. 
Consequently, one can electrically switch the topological edge channels on and off, making the system 
a promising platform for constructing a topological field effect transistor. 


PACS numbers: 73.43.-f, 71.20.-b, 73.20.-r, 73.22.Gk 

I. INTRODUCTION 

The study of topological insulators (TIs) have been one 
of the most active research areas in the past ten years, 
which revolutionized our understanding of the electronic 
band structure. It is now understood that there could 
be nontrivial topologies encoded in the electronic wave- 
functions, characterized by various topological invariants 
according to the symmetry class of the system, and physi¬ 
cally manifested by the appearance of topological bound¬ 
ary states. For example, two-dimensional (2D) TIs, also 
known as the quantum spin Hall (QSH) insulators, are 
characterized by a Z 2 topological invariant and have 
spin helical edge states on sample boundaries,^ for which 
back-scattering is suppressed in the presence of time re¬ 
versal symmetry.Hence they hold great promise for 
applications such as low-dissipation electronics, spintron- 
ics, and quantum computations. For these TIs, the non¬ 
trivial topology as well as the boundary states are pro¬ 
tected by the finite energy gap, i.e., they are robust 
against perturbations as long as the insulating gap does 
not close. 

It was later realized that the topological classification 
could be pushed beyond insulators to states without a 
gap.^“^ In particular, a novel state called Dirac semimetal 
(DSM) has been proposed and successfully demonstrated 
in recent experiments for two crystalline materials NasBi 
and Cd 3 As 2 .^“^^ In these materials, the Fermi energy sits 
at two three-dimensional (3D) Dirac points—where the 
bands touch with a fourfold degeneracy—and the dis¬ 
persion is linear along all three directions in reciprocal 
space. Each Dirac point can be viewed as consisting of 
two Weyl points of opposite chiralities and is protected 
by the crystalline symmetry.Such unusual elec¬ 
tronic structure endows the system with many intriguing 
properties like the surface Fermi arcs and the quantum 
magnetoresistance.Perhaps more importantly, 
DSMs are expected to be an ideal parent compound 
for realizing other novel topological states such as Weyl 
semimetals, TIs, and topological superconductors.^^ Par¬ 


ticularly in this regard, DSMs offer a simple alternative 
to achieve the 2D TI phase through the quantum confine¬ 
ment effect.It has been shown that with increasing 
thickness, DSM thin films exhibit oscillations in the 2D 
Z 2 invariant whenever a quantum well state crosses the 
Dirac point.Hence a QSH phase can be realized by 
a proper control of the film thickness. Since the QSH 
phase has so far been detected in only a few systems, 
given its fundamental and technological importance, the 
new approach to realize it using DSMs would be of great 
interest. Furthermore, the unique properties of DSMs 
may offer new methods to manipulate the properties of 
the QSH phase. 

Motivated by these recent breakthroughs and by the 
great interest in utilizing DSMs for topological devices, in 
this work, we investigate the possibility of electric control 
of the topological phase transitions in a DSM thin film. 
We show that by using a vertical electric field, DSM thin 
films can be switched between a topological QSH phase 
and a trivial insulator phase. Hence one can electrically 
manipulate the topological edge states, and the charge 
and spin conduction through a finite sample can be read¬ 
ily switched on and off. This leads to a simple design of 
a DSM-based topological field effect transistor with ad¬ 
vantages of fast-speed, low power consumption, and low 
dissipation, owing to the robust topological edge channels 
combined with full electric control. 

Our paper is organized as follows. In Sec. II, we dis¬ 
cuss the modeling of a DSM thin film and we present 
an analytic derivation of its low energy levels and its 
topological phase under the influence of a vertical elec¬ 
tric field. In Sec. Ill, we show the results of numerical 
calculation, which explicitly demonstrate the topological 
phase transition and the associated change in topologi¬ 
cal edge states and confirm the physical picture obtained 
from our analytical analysis. Finally, we give a discus¬ 
sion of some aspects of the effect and its possible device 
applications, and summarize our results in Sec. IV. 
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II. MODEL AND ANALYTIC ANALYSIS 

Our analysis is based on a generic low-energy effective 
model describing the DSMs AsBi (A=Na, K, Rb) and 
Cd 3 As 2 as derived in previous works.In these mate¬ 
rials, the states around Fermi energy can be expanded us¬ 
ing a minimal four-orbital basis of | 5 'i/ 2 ,1/2), IP 3 / 2 , 3/2), 
| 5 'i/ 2 , —1/2), and IP 3 / 2 ,—3/2). Around F-point in the 
Brillouin zone, the effective Hamiltonian expanded up to 
quadratic order in the wave-vector k is given by 

Ak^ 0 0 - 

-M{k) 0 0 

0 M{k) -Ak- ’ 

0 -Ak+ -M{k) 

( 1 ) 

where So{k) = Co-\-Cikl-\-C 2 {kl-\-ky), k± = kxhziky, and 
M{k) = -Mo^Mikl^M 2 {kl+kl) with Mo, Mi, M 2 > 0 
to reproduce the band inversion feature at F-point. The 
material-specific parameters A, Ci, and Mi are deter¬ 
mined by fitting the first-principles result or the experi¬ 
mental measurement. It has been shown that this model 
nicely captures the essential low-energy physics as com¬ 
pared with experiment. 

For bulk DSMs, model (1) gives the energy dispersion 
S{k) = eo{k) ± ^M{ky + A/cy, where k\\ = {kx,ky) 

is the 2D wave-vector in the kx-ky plane. The spec¬ 
trum has two Dirac poin ts locate d along the k^-dixis at 
(0,0,±/cd) with /cd = a/Mo/Mi. Each Dirac point is 
four-fold degenerate and can be viewed as consisting of 
two Weyl nodes with opposite chiralities (as represented 
by the two 2x2 diagonal sub-blocks in Hamiltonian 
( 1 )). The dispersion around each Dirac point is lin¬ 
ear in all three directions, as can be seen by expanding 
S{k) at (0,0,r/cD) (r = dz labels the two Dirac points): 

£{k) ^A^k^ -h A^k‘^ -h 4:M^k^{kz — One notes 

that the low-energy spectrum is anisotropic as manifested 
in both the distribution of Dirac points as well as the dif¬ 
ferent Fermi velocities along kz versus that in the kx-ky 
plane (Fermi velocity along kz is typically much slower), 
which leads to quite different behaviors when a DSM is 
confined along different directions. 

DSMs such as Na 3 Bi and Cd 3 As 2 have layered struc¬ 
tures along crystal c-axis. Hence their thin film struc¬ 
tures with confinement along ^-direction can be more 
readily fabricated. Consider a DSM thin film with thick¬ 
ness L confined in the region z G [—L/2, L/2]. For small 
L, the electron motion along ^ will be quantized into dis¬ 
crete quantum well levels due to quantum confinement 
effect. This generally turns the system from a semimetal 
to a semiconductor. Using quantum well approximation, 
each quantum well level has a quantized effective wave- 
vector kz such that {kz)n = ^ and (/c^)n — (riTr/L)^ with 
n{= 1 , 2 , • • •) counting the sub-bands and the angular 
bracket meaning the average over quantum well states. 

One observes that for each sub-band n, Hamiltonian 
( 1 ) has a similar form as the low-energy model describing 


the 2 D QSH systems in HgTe/CdTe quantum wells, 
with the term 

M{k)^ M{n,ki\)=Mn + M2{kl + kl), ( 2 ) 

where Mn = —Mq + Mi{mr/L)'^ is a sub-band depen- 
dent mass which determines the gap of the sub-band at 
F-point of the 2 D Brillouin zone. It’s known that band 
inversion occurs (around k\\ =0) when sgn(A^^/M 2 ) = 
_2 25,26 pg when M-n and M 2 have opposite signs, 
which signals a nontrivial Z 2 = 1 character of the sub¬ 
band n. Given that Mq, Mi, M 2 > 0, this happens when 
kn = (^^)n^^(= nn/L) < /cd is satisfied. Therefore, for a 
thin-enough film such that ki = n jL > k-^, all the sub¬ 
bands are topologically trivial with positive mass terms 
Mn- With increasing film thickness, the system becomes 
nontrivial once the first (n = 1) sub-band has its mass 
Ml inverted when ki < k^- The inverted sub-band con¬ 
tributes a Z 2 = 1 and in the inverted band gap, there ap¬ 
pears a pair of spin-helical edge states protected by time 
reversal symmetry on each edge of the quasi-2D system. 
Further increasing L would invert M 2 oi the second sub¬ 
band, leading to two pairs of edge states. However, for 
Z 2 group: 1 + 1 = 0, hence this state is topologically 
trivial. Physically, it is because backscattering can occur 
between edge states from different time reversal pairs. 
Following this logic, the topological properties as well as 
the bulk band gap show oscillatory behavior as a function 
of the film thickness. 

Since the sub-band dependent mass Mn plays the key 
role in determining the topological properties of the sys¬ 
tem, we shall focus on the change of Mn by a vertical 
electric field, aiming to achieve an electric control of the 
topological phase of DSM thin films. To proceed, one 
notes that the lower diagonal block of Hamiltonian (1) 
is formally the time reversal counterpart of the upper 
block, which share the same energy spectrum and the E 
field does not mix the two.^^ Hence to study the change 
of Mn, it is enough to consider only the upper block 
denoted by h(k). Modeling with the hard-wall bound¬ 
ary condition for the confinement potential, we have for 
sub-band n, 

ft.„(fe||) = £o(n,kii)I+Akxax-Akyay + M(n,kii)az, (3) 

where a’s are the Pauli matrices, I is the 2 x 2 identity 
matrix, and £ 0(^5 ^||) = Cq + Ci(n 7 r/I/)^ + ( 72 ^:^. The 
energy eigenstates are given by 

(r|v|/„„(fe||)) = (4) 

with eigen-energies 

£na(kii) = eo(n,fc||) -hM(n,fe||)2, (5) 

where S is the area of the thin film, a = E, Xna are 
the two eigen-spinors along the quantization direction 


'H(fc) = £o(k) 


M(k) 

Ak_ 

0 

0 
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{Akx,—Aky, M{n,k\\)), and 





nir 

~T 




(6) 


is the quantum well state for the n-th sub-band. The ver¬ 
tical electric field is modeled by adding a diagonal poten¬ 
tial energy term V{z) = eEz where (—e) is the electron 
charge and E is the effective field strength which may be 
considered as including the static screening effects. 

For a qualitative analysis, we assume small field and 
treat V perturbatively. Because V{z) is odd in z, it is 
easy to see that the first order perturbation in energy 
vanishes. The leading order perturbation comes at the 
second order, with 


^ 0.0165), the summation only includes the even integer 
numbers, because V{z) only couples states with opposite 
parities in z. 

Similarly, the energy shift for state |Ti_(0)) can be 
calculated. 


S£i_ 




64L4 1 

7r6 (Ml-Cl)"'’ 


( 10 ) 


which is positive, showing that the coupling induced by 
the E field pushes up the energy of |Ti_(0)). Therefore, 
in the Hilbert sub-space of the first sub-band, the E field 
renormalizes the value of the mass: 


Ml Ml + dMi, ( 11 ) 


(^ 11 ) 


E 




(7) 


where the summation is over all other states different 
from and in reality, it has a physical cutoff for which 
the low-energy description is no longer valid. One notes 
that in order to analyze the renormalized Mn, it is suffi¬ 
cient to focus on the change at T-point of the 2D Brillouin 
zone by setting = 0 . 

We are most interested in the case in which the mass 
(gap) of the first sub-band can be inverted by the elec¬ 
tric field, because then the two sides of the topological 
phase transition have the most salient contrast: absence 
or presence of the topological edge channels, hence lead¬ 
ing to the best on-off ratio when considering a topological 
transistor based on it. For such case, we consider a thick¬ 
ness L such that Mi = —Mq + Mi{7t/L)‘^ > 0, i.e. an 
initially trivial system with Z 2 = 0 in the absence of E 
field. At k\\ = 0, we have Xn+ = | t) Xn- = \ i) foi* 
all n, where | f) and | 1) are the two eigenstates of az. 
For the n = 1 sub-band, we have at k\\ = 0, 


5Si+ « -e^E^ 




m>l 


^m-\- ^ 1 + 


= -e^E^ 


1 


(8) 


where 




m" 


{m? — 1 )^ 


mGeven 


(9) 


In the first equality of ( 8 ) we used the fact that the state 
|Ti+( 0 )) does not mix with the |T^_( 0 )) states from 
the valence bands by the E field because their pesudo- 
spin parts y are orthogonal. Also note that one has 
Afi > I Cl I in order for the model ( 1 ) to describe a 
semimetal phase, hence £m-\- > £i-\- for m > 1 , hence the 
perturbation due to the coupling between |Ti+( 0 )) and 
|T^+( 0 )) (with m > 1 ) generally pushes down the energy 
level of |Ti+(0)), making < 0. In the expression (9) 
of the constant factor 77 (with a rapidly converging value 


with the correction 

^ 2Mi 

~ Ml- 

( 12 ) 

Using this estimation, one observes that the gap of the 
first sub-band would decrease with increasing E field, and 
closes at 


6Mi{E) = 


5£ij^ — 5£i- 


-e^E"^ 


64i 


Ec 


TT^ ImI-cI 

8e2i2 y 27?Mi 


(13) 


which signals a topological phase transition point and 
beyond which the gap reopens with the system turned 
into a QSH insulator phase characterized by Z 2 = 1. 

For large E field with eEL being comparable or even 
larger than Mi^ the perturbative calculation is no longer 
expected to be accurate. Nevertheless, the general phys¬ 
ical pictures from the above discussion still applies: the 
level repulsion due to higher sub-bands would generally 
decrease and invert the gap of the first sub-band, gener¬ 
ating a topological phase transition. We shall explicitly 
demonstrate this in the next section through numerical 
calculations. 

The above analysis can also be applied to higher sub¬ 
bands when a thicker film with the n-th (n > 1 ) sub-band 
most close to transition is considered. We will discuss this 
later in Sec. IV. 


III. NUMERICAL RESULTS 

For numerical investigation, we discretize the model 
(1) on a 3D lattice with lattice constants = 

0.5448 nm (1.264 nm), and with az = 0.4828 nm (2.543 
nm) being set to the interlayer separation pertinent to 
NasBi (Cd 3 As 2 ). The standard substitutions 

1 2 

ki ^ — sm{kiai), k‘f ^ ^[1 — cos(/c^a^)] (14) 

di CLj^ 

are adopted (i = x, y, z) for lattice discretization around 
F-point. Since we require the initial state at = 0 is of a 
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FIG. 1. (color online) Field induced topological phase tran¬ 
sition in NasBi thin film, (a) Variation of energy gap as a 
function of the vertical field E showing the gap closing and 
reopening process, marking a topological phase transition be¬ 
tween Z 2 = 0 and Z 2 = 1 phases, (b-d) Energy spectra 
corresponding to A, B, and C as marked in (a) plotted versus 
kx (with ky — {)). The hrst, second, and third sub-bands are 
marked using red, green, and blue colors respectively. The 
parameters for used in the calculation are I — 5, Co = —63.82 
meV, Cl = 87.536 meV • nm^, C 2 = —84.008 meV • nm^. 
Mo = 86.86 meV, Mi = 106.424 meV • nm^ M 2 = 103.610 
meV -nm^, and A = 245.98 meV -nm. For better comparison, 
a rigid energy shift is applied to make the gap center at zero 
energy. 


trivial insulator phase, we need the number of layers £ < 
[7rMi/(Moa;2)J +1, where [• • • J is the floor function. And 
in order for the band gap to be inverted by a relatively 
small E field, one may wish to have the initial gap size 
2A4i not too large. 

Let’s consider NasBi first. The model parameters we 
use are listed in the caption of Fig. 1, which have been 
extracted from the first-principles calculations and com¬ 
pared well with experiment. For NasBi thin films, the 
critical thickness for which the first sub-band undergoes 
band inversion is around £c — 6.^^ Hence we take a film 
with £ = b layers {L 2.414 nm) for demonstration. 
Fig. 1(a) shows the variation of the band gap Eg as a 
function of the E field. The result is symmetric between 
positive and negative values of E^ so only the positive 
part is shown here. Initially, at = 0, the system has 
a confinement gap about 71 meV (marked by point A). 
With increasing E field, the gap decreases and closes at 
a critical value Ec — 140 mV/nm (marked by point B), 
and then reopens and increases with E. This is consis¬ 
tent with our previous discussion in Sec. II. The value of 



FIG. 2. (color online) Galculated LDOS of a side surface for 
NasBi thin him. (a) is for point A (trivial insulator) and (b) 
is for point G (QSH insulator) as marked in Fig. 1(a). The 
calculation is for a slab which is semi-inhnite along y-direction 
and the parameters are the same as for Fig. 1. 


Ec is also not far from our estimation in Eq.(13) which 
is about 208 mV/nm. We also plot in Fig. l(b-d) the en¬ 
ergy spectrum of the system corresponding to the three 
representative states marked by A, B, and C in Fig. 1(a). 
It shows that on both sides of the gap closing point, the 
system is insulating with the gap belong to the hrst sub¬ 
band. At the critical value Ec^ the band gap closes at 
F-point, marking the topological phase boundary which 
separates the topologically trivial and nontrivial phases. 

To further demonstrate the topological nature of the 
transition and to visualize the edge states, we com¬ 
pute the surface local density of states (LDOS) for the 
side surface. Due to the isotropy in the kx-ky plane of 
the low-energy model (1), without loss of generality, we 
choose the surface perpendicular to ^-direction of the 
quasi-2D system. The surface LDOS p{kx) can be cal¬ 
culated for each kx from the surface Green’s function 
p{kx) = “Tr[ImGoo(^cc)]/7r, where Goo is the retarded 
Green’s function for the surface layer (labled by index 
0) of the lattice.Goo can be evaluated by the transfer 
matrix through a standard numerical iterative method. 
The obtained surface LDOS for states before and after 
the phase transition (for state A and C) are plotted in 
Fig. 2. One observes that for both cases, the confinement- 
induced bulk gap can be clearly identified. Before the 
topological phase transition {E < Ec), there is no states 
inside the gap. In contrast, after transition {E > Ec), 
there appear two bright lines crossing the gap, corre¬ 
sponding to the spin helical edge states for the Z 2 non¬ 
trivial QSH phase. As long as time reversal symmetry is 
preserved, these gapless modes are protected and carriers 
in these channels cannot be backscattered.^’^ Therefore 
transport through these channels is in principle dissipa¬ 
tionless. 

Similar analysis applies to Cd 3 As 2 as well. In Fig. 3, 
we plot the variation of its confinement-induced gap ver¬ 
sus the film thickness, which clearly shows the oscilla¬ 
tion behavior of the gap.^^ One observes that the critical 
thickness £c is at about 37 layers. Here we choose a film 
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FIG. 3. (color online) Confinement induced gap size versus 
thickness i for Cd 3 As 2 thin films. The parameters used in 
the calculation are Co = —219 meV, Ci = —300 meV • nm^, 
C 2 = —160 meV • nm^, Mo = 10 meV, Mi = 9600 meV • nm^, 
M 2 = 180 meV • nm^, and A = 275 meV • nm. 


thickness of ^ = 20 layers {L = 50.86 nm) for demon¬ 
stration. The variations of the gap with respect to the E 
field as well as representative energy spectra are shown 
in Fig. 4. Again the gap closing and reopening process 
similar to Fig. 1(a) is observed. The critical value of 
Ec — 5.26 mV/nm also agrees well with the estimation 
^4.17 mV/nm from Eq.(13). The energy spectra also 
coincide with our previous analysis. Fig. 5 shows the 
side surface LDOS plots for states A and C (marked in 
Fig. 4(a)), clearly showing the appearance of topologi¬ 
cal edge states across the transition. These results show 
qualitatively the same features as those for NasBi. 

Our numerical results discussed above thus confirm our 
analytical analysis in Sec. II. A vertical electric field can 
be used to control the topological phase transitions and 
the topological edge states in a DSM thin film. 


IV. DISCUSSION AND SUMMARY 


This work theoretically demonstrates the possibility to 
electrically control the Z 2 = 0/1 topological phase transi¬ 
tions in a DSM thin film. Since the bulk topology is tied 
to the existence of topological edge channels and hence to 
the charge/spin conductance, this indicates that one can 
achieve a full electric control of the on/off charge/spin 
conductance of such a system, making it a suitable can¬ 
didate for a topological field effect transistor. The electric 
field can be generated by the standard top and bottom 
gates setup. Compared with the traditional MOSFET 
which works by injection and depletion of charge carriers 
in the channel region and has a response timescale de¬ 
pending on factors such as the charge concentration and 
the carrier mobility, the operating mechanism for a topo¬ 
logical transistor is expected to have a high on/off speed 
with electronic response timescale and a better power 
efficiency.In addition, multiple conducting channels in 
a transistor can be achieved by designing a multilayer 
structure with alternating DSM layers and insulating lay- 






- 0.2 0 0.2 
kx/Ti 


FIG. 4. (color online) Field induced topological phase tran¬ 
sition in Gd 3 As 2 thin film, (a) Variation of energy gap as a 
function of the vertical field E. (b-d) Energy spectra corre¬ 
sponding to A, B, and G as marked in (a) plotted versus kx 
(with ky = 0). The first and the second sub-bands are marked 
using red and green colors respectively. The parameters for 
used in the calculation are the same as for Fig. 3 and ^ = 20 
is taken. 



FIG. 5. (color online) Galculated LDOS of a side surface for 
Gd 3 As 2 thin film, (a) is for point A (trivial insulator) and (b) 
is for point G (QSH insulator) as marked in Fig. 4(a). The 
calculation is for a slab which is semi-infinite along ^-direction 
and the parameters are the same as for Fig. 4. 


ers, similar to the structure as in Ref.30. 

For device design, we have seen that a proper film 
thickness can be chosen such that the starting confine¬ 
ment gap is small hence can be easily inverted by a small 
applied field. However, there is a tradeoff because if the 
gap is too small, then the thermally populated carriers 
in the bulk could strongly contribute to the transport. 
Therefore, a balance needs to be achieved for the device 
to have an optimal performance with relatively low power 
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consumption. 

In our analysis, we have focused on the phase transi¬ 
tion in the first sub-band. Similar analysis can also be 
extended to higher sub-bands if a thicker film with its 
n-th (n > 1) sub-band close to phase transition is con¬ 
sidered. For example, consider a film thickness such that 
its second sub-band is just before the gap-closing. In this 
case, we would have AI 2 > 0 and < 0, and the sys¬ 
tem is in a Z 2 = 1 phase. Perturbation to second order 
in the field strength gives the energy correction of 


5 £ 2 + « - e ^£’2 


256i^ 


TfU 


1 


Mi + Ci 

1 


81[3Ci +5Mi -2Mo(L/7r)2] 


, (15) 


for the |^ 2 +( 0 )) state of the second sub-band, with 
v' = Emeodd,m> 3 ^V(^^-4)®. A similar expression 
can be obtained for 882 - as well. The first term in the 
parenthesis of (15) is from the coupling with the m > 2 
sub-bands, while the second term is from the coupling 
with the first sub-band. The sign of this energy shift 
(and hence the correction of M 2 ) would depend on the 
competition between the two terms and is not necessarily 
negative. Nevertheless, for such higher sub-band case. 


even if a field-induced topological phase transition can 
be realized, the topologically trivial phase would in fact 
still possess edge states. Although these (even number 
of pairs of) channels are not topologically robust, their 
presence would make the trivial state not completely ‘off’ 
hence is detrimental to the performance of a transistor. 

In summary, we have demonstrated that full electric 
control of topological phase transitions in a DSM thin 
films can be achieved through the interplay between the 
quantum confinement effect and the coupling between 
sub-bands induced by a vertical electric field. As a result, 
a topological field effect transistor can be constructed 
in which carriers are conducted through the topologi¬ 
cal edge channels. Given that several DSM materials 
have been experimentally demonstrated and that the 
progress in material fabrication technology such as molec¬ 
ular beam epitaxy has allowed film growth with atomic 
precision, it is quite promising for the physical effect and 
the DSM-based topological transistor proposed here to 
be realized in the near future. 
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